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Abstract 

> : 

We perform a systematic study of CP-violating observables in D° /D° decays. 
We show that given a final state, which may be a CP eigenstate or not, all the 
■ observables can in principle be measured if decay-time information is available 

and if, for CP non-eigenstates, both decay rates into / and / are considered. As 
an illustration of our analysis we discuss an effective superweak scenario. 



t^' 1 Introduction 

The study of mixing and decays of neutral D mesons is an important testing ground 
for the Standard Model (SM), with the potential for discovery of new physics. (See, 
for instance, refs. || Q.) In the SM, short-distance contributions to mixing arise 
from the box diagrams [§J, and are negligible ||: Am^ ) ox (S'M) rs 10~ 17 GeV. Initially, 
it was thought that long-distance effects arising from second-order weak interactions 
with mesonic intermediate states would give a much larger contribution, about two 
orders of magnitude larger than the short- distance one ||. It was later argued by 
Georgi that cancellations occur among the various dispersive channels, making 
the overall long-distance contribution smaller than initially estimated. Recently, this 
conjecture has been confirmed by an analysis of D°-D° mixing in the framework of 
heavy quark effective field theory (HQEFT) including leading order QCD corrections 
@, which gave Am" QEFT (SM) « (0.9 - 3.5) • 1CT 17 GeV. One concludes therefore that 
the observation of D°-D° mixing at any existing or planned facility would be evidence 
for physics beyond the SM. 

The study of D°-D° mixing also provides the exciting possibility of detecting new 
sources of CP violation beyond the Kobayashi-Maskawa (KM) mechanism of the SM. 
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Institute for Theoretical Physics, Boltzmanngasse 5, A-1090 Vienna 
t Researcher of the Technical University of Lisbon 
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Indeed, most of the models beyond the SM which predict significant contributions to 
D°-D° mixing contain new sources of CP violation. Recently, it has been emphasized 
that CP violation can play an important role in the study of D°-D° mixing || [H], |ll| . 

In this paper we perform a systematic study of a class of physical observables which 
can provide evidence for CP violation in the D°-D° system. We then assume (as was 
also done in refs. ||, [Tlj ) that the dominant new source of CP violation in the D°- 
D° system arises from an effective superweak interaction (ESI) leading (in a suitable 
phase convention) to a non- vanishing ImM 12 , while and the decay amplitudes are 
dominated by the SM contributions, and are real (we neglect small phases in the KM 
matrix proportional to A 4 in the Wolfenstein parameterization). Although this is a 
specific assumption, it turns out that a large number of models which can lead to 
significant D°-D° mixing belong to the ESI class |I2|. Our analysis is complementary 



to the work of ref. [T3|, where an opposite assumption was made, namely that physics 
beyond the SM contributes significantly only to the decay amplitudes. 

We will start by identifying CP- violating observables in tagged D°/D° decays to 
specific final states, like CP non-eigenstates ///, and CP eigenstates F. We do this 
by studying the behaviour of various quantities under phase redefinitions of the state 
vectors. We point out that the study of the time-dependent decay rates of D° and 
D° to a final state / and its antiparticle state / can, at least in principle, allow the 
measurement of all the CP-violating observables which we have identified. The same 
is true in the case of decays to CP eigenstates. Furthermore, after having switched to 
an ESI scenario, we consider / = K~n + /f = K + n~ and F = K + K~ as illustrative 
examples and possible candidates for an experimental search for CP violation. 



2 Phenomenology of CP- violation in tagged D°/D° 
decays 

We want to describe the decays of D° and D° to a final state \f > and its antiparticle 
state |/ >. We first assume that \f > is not a CP eigenstate. Phenomenologically, 
there are four independent amplitudes 

A(V'^{)=A, A{^^ = X, A(^'^{)=B, A(V'^)=&, (1) 

entering in the description of those decays. The mass eigenstates of the charmed neutral 
mesons are 
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where Am = m H — m L > (i?=heavy, L=light), Ar = T H — T L , and M 12 — |r 12 
is the off-diagonal matrix element in the effective time evolution in the D°-D° space. 
The sign of the square root in eq. @ is chosen such that Am > 0. 

In quantum mechanics, physical quantities cannot depend on the phases of the 
initial and final states. In our case, physics is therefore invariant under the phase 
redefinitions 

\D° >^ ^\D° >, \D°>^e^\D°>, \f >^ e^lf > , |/>-> e ^|/>, (4) 
with four independent phases 7, 7, 7/ and 7/. Noting that under such a redefinition 

A -> e i(7 ~ 7 ' } ,4, 5 -> e'^B, I -> e i( ^ } A, 5 -»■ e'^^B, - -> e i(7 " 7) - , (5) 

we see that the only quantities which are at most quadratic in the amplitudes and 
which have a chance to be physical are 

77, L4I 2 , L4I 2 , \B\ 2 , \B\ 2 , -A*B, V -A*B, -AB*, V -AB\ (6) 

p q p q 

where r\ = \q/p\. Because of the possibility of the phase redefinitions fl5|), there are 
only seven independent quantities in eq. (|6]), five moduli and two phases. One could, 
for instance, eliminate \B\ 2 , \B\ 2 , (q/p)AB* and (p/q)AB* from the set in eq. (|6]) and 
take 77, \A\ 2 , \A\ 2 , \(q/p)A*B\, \(p/q)A*B\, &rg{(q/p)A*B) and arg((p / q) A* B) as a 
set of independent quantities. Note that instead of 77 very often the phase arg(r^ 2 Mi 2 ) 
is used, since 7/ is a function of this phase and of Am, Ar. However, in the following 
discussion this question of independence will not be important. 

In order to study CP violation we first consider the consequences of CP invariance. 
If CP is a good symmetry, there exist two phases (3 and f3f, and a CP eigenvalue (d 
(Cd = sucn that 

cv\v >=] )f} \& >, cr\^' >=]-^\v >, 

cp|{ >=i>*#>, c^>=r^i{>, (?) 

and 

CV\V n > = (d\D h >, 

CV\V C > = -Cd\D l > . (8) 
As a result, the CP invariance conditions are 

A = (SWf-flA, B = e i{ ^ + ^B, 1 = ( D e ip , M 12 = e~ 2ip M* , T 12 = e~ 2ip T* . (9) 

p 

Therefore, if CP is conserved the quantities 
rf - 1, \A\ 2 - \A\ 2 , \B\ 2 - \B\ 2 , X = —A*B - & = —J& - ^AB* (10) 
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all vanish. On the other hand, if CP is not conserved the quantities in eq. (|T0|) measure 
CP violation. Knowing their values means having the complete information about CP 
violation in D°-D° mixing and in the decays of D°/D° into the specific final states ///. 
Furthermore, all observables in this connection must be functions of the quantities in 
the set of eq. (|10D and possibly of 

\A\ 2 + \A\ 2 , \B\ 2 + \B\ 2 , y = —A*B + $ = —M* + j^AB*, (11) 

V V 

which are not CP-violating. CP-violating observables have to be odd functions of the 
quantities in eq. (|T0l). 

It is easy to see that X is related to X and y by 

X = -\{v 2 + \w + ^-- e )T. (12) 
2 r] z G r/ fc 

As we will see below, instead of using Re^ and Im^ (or KeX and ImA?) as independent 
CP-violating observables, it is more convenient to take 

Im [(Am - -Ar)*] = -2 Im [M 12 (AB* + A*B)) - Re [T 12 (AB* - A*B)] , 

Im [(Am - ^Ar)^] = -2 Im [M 12 (AB* + A*B)} + Re [T l2 (AB* - A*B)} . (13) 



It can be checked by using eq. ( [12]) that both sets are equivalent if AmAr ^ 0. To 
derive eqs. flTBp we have used eq. (^). 

It is important to realize that there are three possible different kinds of CP violation 
in the D°-D° system (or any other neutral- meson system). Namely, there may be CP 
violation in mixing (which just means i] ^ 1), direct CP violation (that is, \A\ ^ \A\ 
or \B\ 7^ \B\), or CP violation in the relationship between mixing and decays — non- 
vanishing X or X. Often, authors refer to this last form of CP violation as "CP 
violation in mixing", but we think this designation is inappropriate. CP is violated 
in mixing if and only if rj is different from 1, and this is independent of the decay 
amplitudes to any particular decay modes. 

If at t — a D° or a D° is produced, then its time evolution is given by 

\D°(t)> = f + (t)\D >+«f_(t)\D >, 

V 

\D°(t)> = ^f4t)\D°>+f + (t)\D°>, (14) 

q 

m = ^(m H + m L ) , T = ^{T H + T L ) . (16) 

These formulae allow us to consider the time-dependent decay rates of D°/D°. Now 
we address ourselves to the question of whether one can measure, at least in principle, 



with 



where 
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all the CP- violating observables in the set (|10D , which we derived from rephasing in- 
variance, by considering the above-mentioned time-dependent decay rates. The answer 
is affirmative. In order to see this, consider the quantities IT( 



r(D u (i) -> f) - r(D u (t) - /) = e- rt {\A\* -\A?+ (H) Im [(x - iy)X}+ 



+i(rt) 2 \(v 2 \b\ 



\B\ 2 )(x 2 + y 2 ) - {\A\ 2 - \A\ 2 ){x 2 - y 2 ) + £>(U 3 )} 



T(D°(t) -> /) - T(D°(t) -> /) = e- rt {\B\ 2 



\B\ 2 + (Tt) lm[(x-iy)X] + 



+ \ (Tt) 2 U V 2 \A\ 2 - MA\ 2 )(x 2 + y 2 ) -{\B\ 2 - \B\ 2 ){x 2 



0(U*)}, (17) 



with 



x 



Am 



Ar 
2F 



The terms with time dependence e~ Tt are \A\ 2 — \A\ 2 and \B\ 2 — |-B| 2 , whereas those 
with time dependence (Tte~ rt ) coincide with the quantities in eq. (13). 

Taking also into account the sum of the time-dependent decay rates, the moduli of 
all the amplitudes are found from the terms with time dependence e~ rt . This would in 
principle allow the determination of rf from the {(Tt) 2 e~ rt ) contributions in eqs. (|17D. 
However, this contribution can rather be used to get information on x and y, while 77 
is better extracted from inclusive semileptonic decays of tagged D° / D° by 



M--M+ _ T] 4 -l Im (r* 2 M 12 ) 
A/l+A/; ~ ?? 4 + l ~ |M 12 | 2 + i|r!, 



(19) 



where A/+ (A/"_) is the number of positively (negatively) charged leptons originating 
from D°(t) (D°(t)). Of course, with like-sign dilepton events from correlated D°D° 
decays (r/ 4 — l)/(^ 4 + 1) can also be measured. 

If we consider the decays to CP eigenstates \F > the analysis becomes simpler. We 
must then make the identification 



B —> A, B —> A 
in the previous formulas. Then we have X = & and 

Im [(Am - ^Ar)^] = - A Im (M^Xfc) . 



(20) 
(21) 

We note that eq. ( pT|) does not depend on Ti 2 and the two eqs. (|T7| ) become identical in 
the case D°/D° — ► F . Instead of the quantity on the right-hand side of eq. (21) we could 
take Im (^A*A) as the CP- violating quantity which is zero if CP is conserved. Indeed, 
if CP is conserved, there exists a CP eigenvalue £f 
with 

A = (f e~ if5 A . 



±1 such that CV\T >= (^J 7 >, 

(22) 



Therefore, 



are both real. 



q 
v 



A* A = CdCfI^I 2 > M 12 AA* = (M 12 AA* 



(23) 
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3 The effective superweak interaction 



In the SM, CP violation in the decays of D°/D° is extremely small. This can easily be 
seen when using the Wolfenstein parameterization of the KM matrix. Going beyond 
the SM, we assume that, either we can neglect the new contributions to the decay 
amplitudes and to r 12 , or the new contributions to the decay of the c quark are real 
relative to the SM ones. On the other hand, the additional interactions enhance M 12 = 
|M 12 |e 4<?i . Then, if the phase is large, we have an effective superweak scenario with 
| M12 1 3> | r i2 1 , while II [Tl| all the CP- violating quantities in eq. ([U]) arise only from 
0. Assuming for simplicity that Ti2 is exactly real and that there is no CP violation 
in A, B, A, and we obtain 

Am«2|Mi 2 |, Ar^2r 12 cos0, (24) 

and 

\A\ 2 - \A\ 2 = \B\ 2 - \B\ 2 = 0, 
Im [(Am - ^Ar)*] « Im [(•$ - 1— )&] « -G Im [M^AB* + A*B)} , 

« 2^ sin 0. (25) 
7/* + 1 Am 

Since in an ESI the phase is the only source of CP violation, there must be relations 
among the quantities of eq. (fT0|), as borne out by eq. fl25|). 

Considering first CP violation in mixing, we must take into account that eq. (|P9|) 
is somewhat misleading, because it hides the fact that A+ and M- are suppressed by 
mixing, and really contain a factor x 2 + y 2 ~ x 2 . Taking this into account, we see that 
in order to detect CP violation in mixing one has to overcome a suppression factor 
AmAr/r 2 . In the SM we expect AT/T ~ Am{SM)/Y ~ 10~ 5±1 ; therefore, even with 
x 2 ~ 10~ 2 at the present upper limit and sin ~ 1, it seems to be hard to see an effect. 

Taking this crude argument into account, we observe that with an ESI only the 
terms linear in t inside the braces of eq. ([H]) are relevant: 

rp°(f) - /) - r(D°(t) -> /) « r(2?°(t) - /) - r(5°(t) - /) 

» -2(rte~ r ') |AB|xcos5sin0. (26) 

The phase 6 = arg(AB*) is induced by different final-state interactions in the ampli- 
tudes A and B, and is estimated to be ITTH between 5° and 13°. Therefore, cos 5 ~ 1. 



In the case of CP eigenstates F, there is no such phase and 

T{D°{t) -> F) - T{D°(t) ^ A) ^ -2 {Ft e - rt ) \A\ 2 ( F x sin . (27) 

The best final-state candidates to search for mixing and CP violation are probably 
/ = K~n + /f = A+vr", with branching ratios BR(£>° -> A^tt+) = (4.01 ± 0.14) • 
10~ 2 and BR(£>° -> A+vr") = (3.1 ± 1.4) • 10~ 4 , respectively^ and F = K + K~, with 

lr This amounts to CP conservation in the interactions relevant for the decays, with a phase con- 
vention (3 = f3f = in eq. (Q). 

2 See, however, the remarks to the latter value in refs. [0, n4|. 
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Cf = 1 and BR(D° -> K+K~) = (4.54 ± 0.29) • 10~ 3 . The sensitivity of these final 
states to CP violation is nearly the same, because 

mm 



BR(K-tt+)BR(K+7v- 



For f — K iT + /f = K + tt the decay rates are usually normalized to the number of 
Cabibbo-allowed decays. We thus consider the distribution in decay time 

[v(D°(t) - /) - r(D°(t) -> f)]dt ^ [r(D°(t) -> /) - rp°(0 - ^ 



f °°[T(D%t>) - /) + rpo(t') - jnr(L>°(f) - /) + r( J oo(t') - /)]<# 

r D o - k-J) x cos 5 sin * rt e ' ( 29 ) 

Experimentally, the bound on x is roughly 1/10. Therefore, with 10 3 — 10 4 Cabibbo- 
allowed decays, or 10-100 doubly- Cabibbo-suppressed decays, and with decay-time 
information, it might be possible to constrain xsin0 beyond the trivial constraint 
|xsin0| < x. The effect of CP violation with an ESI in the determination of mixing 
from T(D°(t) -»• f) + T(D°(t) -> /) has been studied in refs. [0, |TJ. It was found 



that, at the present level of accuracy, it is included in the experimental errors |il 
For F = K + K~ we use the normalization 

[T(D°(t) —> F) — T(D°(t) — > F)\dt ^ 



f™[T(D°(t') -> F) + T(D°(t') -> 



xsin0rte _lt d(rt). (30) 



Here we estimate that 10 2 — 10 3 K + K events would probably give a useful constraint 
on xsmd). 



4 Conclusions 

We have done a systematic study of CP violation in tagged D°/D° decays, in which 
we have identified all the rephasing-invariant observables which measure CP violation 
in the decays to a final state / and its antiparticle state /, or to a CP eigenstate 
F. These quantities, which parameterize the complete information on CP violation 
in this context, are constructed from the mixing parameter q/p and from the decay 
amplitudes, and therefore depend on the specific final state. We have shown that by 
studying the time- dependent decay rates of D°/D° to /// or to F one can, in principle, 
measure all the CP-violating observables which we have identified. The D°-D° system 
is unique (see 0] for a nice presentation of this fact), in comparison to the K°-K° and 
B°-B° systems, in that x and y are both much smaller than one. As a consequence, 
one can expand the time-dependent decay rates in terms of (Tt) n e~ rt (n = 0, 1, . . .), 
and it is probably sufficient to keep only terms with n < 2. It is crucial for our analysis 
that experimentally the coefficients of the terms with n < 2 can be extracted. 

Although our analysis is model-independent, we have considered the special case 
of the ESI hypothesis, which assumes that physics beyond the SM does not contribute 
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significantly to the D°/D° decay amplitudes, but gives the dominant contribution to 
M 12 , with an unsuppressed CP- violating phase = argM 12 . Under this assumption, 
we have tried to assess the sensitivity of experiments with decay-time information to 
CP violation. Thereby we have compared CP non-eigenstates (where one has Cabibbo 
non-suppressed versus doubly-Cabibbo-suppressed amplitudes) with CP eigenstates 
(single Cabibbo suppression). We have seen that experimental information on CP 
violation restricts (or determines) the quantity || x sin 0, and the upper limit obtained 
for this quantity must be smaller than the upper limit on x, in order to get a non- 
trivial restriction in the x — <fi plane. In the case of x close to its upper bound, we 
have estimated the numbers of events necessary for this purpose. They seem to be 
well within the range of current and future search for D°-D° mixing ||. Of course, 
an experimental determination of x sin <fi would not only prove the existence of CP 
violation in the D°-D° system, but also give a lower bound on the mixing parameter 
x. 
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